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ABSTRACT 

We demonstrate that rapid expansion of the shocked surface layers of an O- 
Ne-Mg core following its collapse can result in r-process nucleosynthesis. As 
the supernova shock accelerates through these layers, it makes them expand so 
rapidly that free nucleons remain in disequilibrium with a-particles throughout 
most of the expansion. This allows heavy r-process isotopes including the ac- 
tinides to form in spite of the very low initial neutron excess of the matter. We 
estimate that yields of heavy r-process nuclei from this site may be sufficient to 
explain the Galactic inventory of these isotopes. 



Subject headings: nuclear reactions, nucleosynthesis, abundances 
— supernovae: general 



shock waves 



1. Introduction 

In this Letter we propose a new astrophysical site for r-process nucleosynthesis. In 
particular, we show that heavy nuclei with mass numbers A > 130 through the actinides 
may be produced in the shocked surface layers of an O-Ne-Mg core following its gravitational 
collapse. Stellar evolution models show that stars of ~ 8-1 IMq develop degenerate O- 
Ne-Mg cores, at least some of which eventually collapse to produce supernovae (SNe; e.g.. 



Nomoto 1984, 1987: Garcia-Berro fc Iben 1994; Ritossa et al. 1996, 19991: Garcia-Berro et al. 
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19971 : llben et al.l 119971 : iPoelarends et al.l 120071 ). Such a core is separated from the hydrogen 



envelope by a thin C-0 shell and an even thinner He shell. However, the density falls off 
so steeply in the region immediately above the core that we shall refer to the C-0 and He 
shells as the surface layers of the core. These layers result from He and H burning coupled 
with convection during the pre-SN evolution, and this pre-supernova burning can give rise 
to a metallicity-independent neutron excess through reaction sequences such as ^^C(p, 7)^^N 
followed by (3^ decay of ^^N to ^^C (Nomoto 2007, personal communication). The SN shock 
rapidly accelerates as it propagates through the surface layers of the core. This gives rise to 
fast expansion of the shocked ejecta on timescales of ~ 10^^ s. Together with an entropy 
of 5 ~ 100 (in units of Boltzmann's constant k per nucleon) and an initial electron fraction 
of Ye ~ 0.495 (e.g., for a composition of ^'^C : ^^C : ^^O ~ 1 : 3 : 3 by mass), this fast 
expansion enables an r-process to occur in the shocked ejecta, producing neutron-rich nuclei 
with A > 130 through the actinides. 

We discuss the conditions during the expansion of the shocked ejecta in ^ Nucleosyn- 
thesis calculations for two different initial compositions are presented in ^with and without 
neutrino reactions on free nucleons and a-particles. We discuss the implications of this 
new r-process model for abundances in metal-poor stars and for general Galactic chemical 
evolution in ® 



2. Conditions During Expansion of Shocked Surface Layers 



In the generic scenario of stellar core collapse, the inner core becomes a proto- neutron 
star (PNS) and bounces on reaching supra-nuclear density. This launches a shock, which 
propagates into the still collapsing oute r core and falters. According to the neutrino-driven 
SN mechanism (IBethe fc Wilson! Il985l ). the shock is re-energized through heating of the 
material behind it by the neutrin os emitted from t he PN S. Using Nomoto' s (1984, 1987) 



model for a 1.38 M(7) 0-Ne-Mg core , iMayle fc WilsonI (119881 ) , and most recently, iKitaura et al. 



(120061 ) (see also lJanka et al.ll2007l ) indeed obtained neutrino-driven explosion although with 
very different final explosion energies. For our purpose, the exact explosion energy does not 
matter so long as the re-energized shock propagates through the surface layers of the core 
with a sufficient speed. As only stars in a very narrow range of ~ 8-11 Mq develop 0-Ne-Mg 
cores, we consider Nomoto's (1984, 1987) 1.38 core model, which was evolved from the 
He core of a ~ 9 Mq progenitor, to be representative and adop t its quantitative descrip tion 
of the pre-SN c ondit ions. In addition, we use the simulations of iMayle &: WilsonI (119881 ) and 
Kitaura et al.l (120061 ) based on this model as a guide to the shock propagation. 



At the onset of collapse, the core is surrounded by ~ 0.1 Mq of C and O, which in turn is 
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surrounded by several 10~^ Mq of He. We are interested in the outer C-0 layer that contains 
several 10~^ Mq. This layer is located at radii r ^ 10*^ cm and has a thickness of ~ 5 x 10^ cm. 
The density p within this region changes from 10^ to 10^ g cm~^. The steepness of this fall- 
off is measured hy w = — (ilnp/dlnr, which is 61.7 at p = 3 x 10^ g cm~^. This can be 
understood from the condition of hydrostatic equilibrium, —{l/p)dP/dr = GMc/r^ , where G 
is the gravitational constant, Mc = 1.38 Mq is the core mass, and P is the pressure provided 
by electrons and nuclei with temperature T ~ 6 x 10^ K. The pre-SN density structure 
changes as a result of core collapse. However, the local free-fall timescale at r 10® cm is 
~ a/ r^/ (GMc) ~ 0.1 s, comparable to the time between the onset of core collapse and shock 
arrival at this radius. To zeroth order, we assume that the pre-SN density structure in the 
region of interest remains unchanged before being shocked. 

The SN simulations mentioned above showed that the shock speed in the region of in- 
terest is Vsh ~ 10^° cm s^^. For such Vsh, the strong shock condition pw^j^ ^ P is satisfied 
and the energy density of the shocked matter is expected to be dominated by the contri- 
butions from relativistic particles (radiation and electron-positron pairs). In this case, the 
density, velocity, and pressure o f the shocked rn atter (with the subscript "p" standing for 
"post-shock" ) are given by re.g.. lchevali^ll976l ) 



Pp = 7p, Vp = ^Wsh, Pp = ^pt'sh. (1) 

respectively. Using Eq. ([T]) and Pp = (ll/12)ai.ad^p with a^ad being the radiation constant, 
we obtain the post-shock temperature 

Tp = 1.05 X W'p'/'vZo K (2) 
and the post-shock entropy in relativistic particles 

^ = ll'^I^-P. = 56.1^^ k nucleon-\ (3) 
■J ^^aPp pg' 

where pe is the pre-shock density in units of 10^ g cm~^, fsh,io is the shock speed in units 
of 10^° cm s~^, Na is Avogadro's number. The entropy S is also a measure of the energy 
density in relativistic particles relative to that in nonrelativistic ones. Equation gives 
S* 3> 1 in the region of interest, so the corresponding energy density of the shocked matter 
is indeed dominated by relativistic particles. 



Matzner &: McKed ( 1l999l ) showed that Vg^ at a specific r depends on the pre-shock den- 
sity at r and the amount of matter entrained by the shock prior to reaching r. As the region 
of interest is very thin and adds very little to the matter already entrained by the shock, we 
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expec t that the shock speed in this region takes the form Vsh oc p (iMatzner &: McKee 



19991 ). which gives 



^ = -0-19f^l^ = 0.19«;f^). (4) 

dr \dmr J r \ r J 

We assume that matter moves with a constant velocity Vp after being shocked. Consider a 
mass element initially sandwiched between radii r and r + Sr. When the matter at r + 5r is 
shocked, the matter initially at r has moved to r + {vp/vsh)Sr = r + (6/7)5r. The reduction 
of the mass element's thickness from 5r to 5r/l results in the factor of 7 increase in its 
post-shock density (see Eq. [1]). Subsequently, the matter at the inner and outer radii of the 
mass element move with velocities of Vp and Vp + {dvp/dr)6r, respectively. It can be shown 
that the density p of the shocked mass element evolves with time t as: 

^^'^ ^ [i + it/nW[i + it/r,)y 

where t = is the time at which the mass element is shocked, and 

Tl = — T2= ^ = (6) 

Vp' 7{dvp/dr) 1.33w 

In the denominator of Eq. ([5]), the first factor accounts for the expansion of the surface area 
of the mass element and the second for the increase of its thickness. In relating T2 to ti in 
Eq. ([6]), we have used Eqs. ([1]) and (jlj). To very good approximation, the expansion of the 
shocked mass element is adiabatic. This then specifies the time evolution of its temperature 
T as 

T 

T(f) = (7) 

^' [l + (t/ri)]2/3[l + (t/r2)]i/3- ^' 

Considering the proximity of the region of interest to the PNS, we also study how 
interaction with the intense neutrino flux affects nucleosynthesis. Due to the short timescales 
Tl ~ 10^^ s and T2 ~ 10^^ s for the expansion of a shocked mass element, the neutrino 
luminosity and energy spectra do not change during the expansion. We take the luminosity 
to be 4 X 10^^ erg s~^ per neutrino species, {El) / {E^,) = 11 and 16 MeV for z/g and Ue, 
respectively, and {E^) = 25 MeV for = v^, z/^, z/,-, where {E^) and {E^) are the 
averages of the neutrino energy and its square, respectively, over the relevant spectrum 
(the rates of interest for z/g and z/g are essentially determined by the values of {El) / {E^) 
for these species and insensitive to their detailed spectra, while the rates of interest for 
assume a black-body-like spectrum with a temperature of 8 MeV). We include the reactions 
Pe + n ^ p + e~, Ue + P ^ n + e+, + "^He — > ^He -|- n, and + ^He — > -|- the rates 
of which are 0.253, 0.250, 1.28 x 10"^, and 1.40 x 10"^ s"^ per target, respectively, at a 
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radius of 10^ cm (iQian fc Woosleylll996l : IWoosley et al.lll990l ). For a mass element initially 
at radius r, its radius f evolves as 



f(t) =r[l + (t/ri)]. 
All the neutrino reaction rates for this mass element scale as [f(t)]~^. 



(8) 



3. Nucleosynthesis in Shocked Ejecta 

We focus on a specific mass element initially at r = 10^ cm. The pre-shock density 
structure at this position is characterized by pg = 0.3 and w = 61.7. We take fsh,io = 1-5 at 
this position, which gives Tp = 9.52 x 10^ K, n = 7.8 x 10"^ s, and rs = 9.48 x 10"^ s. The 
postshock density is pp = 2.1 x 10^ g cm~'^, corresponding to an entropy of S' = 139. We 
study the nucleosynthesis during the expansion of the above mass element for two different 
initial compositions: (1) X(^^C) = X(^^O) = 0.5 corresponding to Vg = 0.5 as in Nomoto's 
(1984, 1987) original model, and (2) Xi^^C) = 0.14, X{^^C) = X{^^0) = 0.43 corresponding 
to Ye = 0.495, which allows for possible production of ^^C due to mixing during the pre-SN 
evolution. Here, for example, X(^^C) is the mass fraction of ^^C. The above two cases are 
referred to as trajectories 1 and 2 when neutrino reactions are not included and as trajectories 
1' and 2' when these reactions are included. 

To study the nucleos ynthesis associated wit h trajectories 1 and 2, we use the Clemson 



nuclear reaction network (j Jordan fc Meyen l2004j ) . which includes nuclear species from neu 



trons and protons through the actinides. Because the network flow reached the top of the 
reaction network in some cases, we included nuclear fission by allowing the most massive 
species in the network (^'^^U) to fission into two fragments (with Z = 42 and 50, respec- 
tively). This is a crude treatment that will require improvement in follow-up calculations. 
Nevertheless, it illustrates the effect of fission cycling. 

Figure [T] shows the final abundances for trajectories 1 (top panel) and 2 (bottom panel) 
as a function of A. Trajectory 1 (initial Y^. = 0.5) makes heavy, neutron-rich nuclei, although 
the abundance pattern does not resemble the solar r-process abundance distribution well. 
Trajectory 2 (with initial Yg = 0.495) makes heavy, neutron-rich nuclei, including actinides, 
and the resulting abundance pattern shows three distinct peaks in rough agreement with the 
solar r-process distribution. The peaks made by trajectory 2 lie several mass units above the 
corresponding solar r-process abundance peaks. This is mainly due to the rather extensive 
exposure to neutron capture for this trajectory: there are still about 15 free neutrons per 
heavy nucleus at the end of the run. The reason why trajectories 1 and 2 are able to produce 
heavy, neutron-rich nuclei despite the zero to low initial neutron excess is the persistent 
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nucleon-a-particle disequilibrium discussed by lMeyerl ( 120021 ) . In these expansions, the initial 



distribution of nuclei is quickly broken down into nuclear statistical equilibrium, which is 
dominated by free neutrons and protons. As the material expands and cools, the free nucleons 
assemble into a-particles and heavier nuclei. Because of the rapidness of the expansion, 
however, the free nucleons do not assemble into heavier species (particularly a-particles) as 
quickly as equilibrium demands. This leaves a large excess of these free nucleons, which push 
the heavy nuclei to very high mass. For example, for trajectory 2, at Tg ^ 4, the abundance 
distribution is dominated by free neutrons and protons, a-particles, and nuclei with mass 
number A ^ 140. Specifically, the mass fractions of free neutrons and protons, a-particles, 
and heavy nuclei at this temperature are ~ 1.3%, 0.3%, 97.6%, and 0.6%, respectively. The 
heavy nuclei serve as the seeds for the subsequent r-process nucleosynthesis. 

We expect that the final abundance distribution will depend on the details of the tra- 
jectories of the mass elements at late times. Indeed, we note that for trajectory 2, expansion 
occurs so rapidly that reaction freezeout occurs because the density drops to such a low 
value that neutron-capture timescales become long, not because the neutrons are all con- 
sumed. The result is that there are still about 15 free neutrons per heavy nucleus at the 
end of the run (the final mass fractions of free neutrons, a-particles, and heavy nuclei are 
~ 0.2%, 98%, and 2%, respectively). We confirmed this by running a calculation identical 
to trajectory 2 except that the expansion slowed at late times and allowed the neutrons to 
be all consumed. The resulting pattern was broadly similar to that in the lower panel of 
Fig. [H but the contrast between the peaks and valleys was smaller. 

We also followed trajectories 1' and 2', which were identical to trajectories 1 and 2 except 
that neutrino reactions were included. The neutrino interactions had negligible effects on 
the abundance patterns. This is because expansion carried the material away so rapidly 
that the total number of neutrino interactions per nucleon or nucleus was <^ 1 during the 
expansion. 



4. Discussion and Conclusions 



We have proposed a new model of r-process nucleosynthesis in 0-Ne-Mg-core col- 
lapse SNe. Unlike previous models base d on assumed extremely neutron- rich ejecta (e.g.. 



Wheeler et al. 1998 



Wanajo et al.l l2003l ). our model relies on the fast expansion of the 



shocked matter in the weakly neutron-rich surface layers of an 0-Ne-Mg core to achieve 
a dise quilibrium between free nucleons, a-particles, and heavier nuclei. As shown by iMeyer 
(I2OO2I ). the significant presence of free nucleons in disequilibrium facilitates the production 
of seed nuclei with A ~ 140 during the expansion. For matter with an initial Yg = 0.495, the 
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neutron excess is sufficient for an r-process with fission cycling to occur, producing domi- 
nantly the nuclei with A > 130 through the actinides. The key to the fast expansion of the 
shocked matter is the steep density gradient in the surface layers of an 0-Ne-Mg core, which 
enables the inner and outer surfaces of a mass element to expand with significantly different 
velocities, thus making its density drop much faster. 

0-Ne-Mg-core collapse SNe were also proposed as the source for the heavy r-process 
elements (r-elements) with A > 130 (Ba and higher atomic numbers) based on considera - 



tions of Galactic chemical evolution (e.g., iMathews et al.lll992l : llshimaru fc Wanaidll999l ). 



A strong argument for this prop osal was presented using observations of metal-poor stars 
( IQian &: Wasserburgl |2002| . |2003| ). In addition, the key issues regarding the sources for the 
r- elements based on ob s ervat ions and basic understanding of stellar models were sumraarize d 
in bian fc Wasserburgl (|2007|). Data on the metal -poor stars CS 31082-001 JhIU et al.lbooi ). 



HD 115444, and HD 122563 (IWestin et all 120001 1 show that their abundances of the heavy 
r-elements differ by a factor up to ~ 10^. In contrast, these stars have essentially the same 
abundances of the elements between O and Ge (e.g., [ Fe/H] = log(Fe/H)— log(Fe/H)f^ ~ —3). 
Furthermore, when CS 22892 -052 ([Fe/H] = — 3.1, ISneden et al.l l2003l ) is compared with 
HD 221170 ([Fe/H] = -2.2, llyans et all \20od] a nd CS 31082-001 ([Fe/H] = -2.9) with 
BD -M7°3248 ([Fe/H] = -2.1, ICowan et al1l2002h . data show that the stars in either pair 



have nearly the same abundances of heavy r-elements but the abundances of the elements 
between O and Ge differ by a factor of ~ 8 and 6 for the former and latter pair, respectively. 
These results appear to require that the produ ction of the heavy r-elements be deco upled 
from that of the elements between O and Ge (IQian &: Wasserburgl |2002| . |2003| . 120071 ). As 
Fe-core collapse SNe from progenitors of > 11 Mq are the major source for the latter group 
of elements at low metallicities, this strongly suggests that such SNe are not the source 
for the heavy r-elements. The elements between O and Ge are produced between the core 
and the H envelope by explosive burning during a core-collapse SN or by hydrostatic burn- 
ing during its pre-SN evolution. Models of 0-Ne-Mg-core collapse SNe show that the total 
amount of material ejected from between th e core and the H envelope is only ~ 0.01-0.04 
( iMayle &: Wilsonlll988l : iKitaura et al.ll2006l ). much smaller than the ~ 1 Mq for Fe-core col- 
lapse SNe. Thus, 0-Ne-Mg-core collapse SNe contribute very little to the elements between 
O and Ge. The decoupling between these elements and the heavy r-e lements can then be 
explained by attri buting the heavy r-elements to such SNe as argued in lOian fc Wasserburg 
( I2OO2I . I2OO3I . I2OO7I ). This attribution is supported by the r-process model proposed here. 



The collapse of 0-Ne-Mg cores may comprise up to ~ 20% of a l l core -collapse events 
and may thus occur roughly once a century (e.g., iPoelarends et al.l 120071 ). Assume that 
this frequency of occurrence is associated with a reference gas of ~ 10^*^ Mq and that the 
0-Ne-Mg-core collapse SNe over the Galactic history of ~ 10^° yr enriched this gas with a 
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solar mass fraction of ~ 4 x 10~^ of heavy r-elements. We find that each SN must produce 

~ 4 X 10~^ Mq of heavy r-elements. The total mass fraction of heavy r-elements produced 
during the expansion associated with trajectory 2 is ~ 0.02 and a total ~ 3 x 10~^ Mq of 
ejecta may experience such expansion in Nomoto's (1984, 1987) 0-Nc-Mg core model. This 
assumes an initial composition with Yg = 0.495 and would give ~ 6 x 10^^ Mq of heavy 
r-elements per SN. A slightly smaller Ye = 0.49 would double the initial neutron excess and 
the final r-process production, giving ~ 10~^ Mq of heavy r-elements. Thus, we consider 
the model proposed here has the potential to account for the Galactic inventory of such 
elements. To further test this model, we urge two lines of important studies: (1) calculating 
the evolution of ~ 8-11 Mq stars to determine the pre-SN conditions of 0-Ne-Mg cores, 
especially the neutron excess and density structure of the surface layers; and (2) simulating 
the collapse of such cores and the subsequent shock propagation to determine the conditions 
of the shocked surface layers. As these layers contain very little mass, simulations with 
extremely fine mass resolutions arc required to demonstrate the fast expansion of shocked 
ejecta that is proposed here as the key to the production of heavy r-elements. The pre-SN 
abundance of ^^C in the surface layers also requires investigation as this appears critical to 
the total yields of heavy r-elements. 

We thank Ken Nomoto for providing his 0-Ne-Mg core model and for valuable discus- 
sion. We also thank Jerry Wasserburg and an anonymous referee for helpful comments. This 
work was supported in part by DOE grant DE-FG02-87ER40328 and NASA Cosmochemistry 
Program grant NNX07AJ04G. 
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Fig. 1. — Final abundances versus mass number for trajecto ries 1 (top panel) and 2 (bottom 
panel). The (arbitrarily scaled) solar r-process abundances ( iKappeler et al.lll989l ) are shown 
as x's for comparison. The final mass fractions resulting from both trajectories are ~ 98% 
a-particles and ~ 2% heavy nuclei. 



